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Abstract — In this paper, bit error rate (BER) expressions 
for the i-MRC receiver for a decode and forward (DF) 
cooperative system are obtained for Nakagami-m fading, 
where m is not an integer. Previous results were available 
only for integer values of m. BER analysis is done by 
employing approximate statistics of a gamma conditionally 
gaussian (CG) random variable (RV) obtained through the 
Loskot-Prony approximation. Numerical results obtained 
using the analytical BER expressions are shown to closely 
follow the simulation results, despite the cumulative distri- 
bution function (CDF) of the gamma CGRV being a high 
signal to noise ratio (SNR) appoximation. 

Index Terms — BER, Gamma CG distribution, 
Nakagami-m fading 
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Fig. 1. Three node cooperative diversity system. 



I. Introduction 

BER expressions for DF cooperative systems, 
compared to amplify and forward (AF) cooperation, 
are difficult to evaluate, and hence there is con- 
siderable interest in finding analytical expressions 
for the BER for DF cooperative systems. Exact ex- 
pressions for the BER for the piecewise linear (PL) 
combiner were first obtained for Rayleigh fading 
in [1] and [2] for noncoherent binary frequency 
shift keying (BFSK) and binary phase shift keying 
(BPSK) respectively. Results for the more general 
Nakagami-m fading were first obtained in [3] using 
the approach in [2] followed by a simpler approach 
in [4]. [4] also included BER expressions for the 
/l-MRC receiver proposed in [5]. 

One common feature of the above literature is 
the restriction of the Nakagami fading paramter m, 
to being an integer. To the best of our knowledge, 
there has not been any attempt to evalute the BER 



for popular DF receivers like the /l-MRC or PL 
combiner for noninteger values of m. 

In this paper, we use a Loskot-Prony [6] approx- 
imation for the CDF of a gamma CGRV to obtain 
the BER for a /l-MRC cooperative system. This 
approximation for the CDF is known to be tight 
for high SNR. However, through numerical results, 
we show that the related expression for the BER 
obtained in this paper using this approxmiate CDF, 
exactly follows the simulation results. 

In the beginning section, the system model is pre- 
sented, followed by BER analysis. Numerical and 
simulation results are discussed next. Our conclu- 
sions are summarized in the final section, outlining 
the scope for future work. 

II. System Model 

The classic three node cooperative system in 
Figure 1 is considered, where, without loss of gen- 



erality, h represents the Nakagami-m channel gain 
with fading figures m and Q., E the transmit power at 
a node, x the transmitted symbol at a node, and the 
subscripts s and r the source and relay parameters 
respectively. 

A. A-MRC 

The decision statistic for the i-MRC receiver for 
BPSK modulation, is given by [4], [5] 



X + AY ^ 0, 0<A<1 
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where X ~ N (a s h 2 s ,b s hf), Y ~ N(a r h 2 r ,b r h 2 ) with 

a > ~ iff' ^' = ~No' C > = Hi ' ' G ^' ^' 

p h 2(x) = — — exp (-Cjx) , x, Ci > 0, m, > 0.5 
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^ ~ G(ci,rtii), where ^ denotes the Gamma dis- 
tribution [8]. Assuming equal probability of the 
transmitted symbol x s = {1,-1}, from (1), the 
average BER for a /1-MRC cooperative system can 
be expressed as 



xMl-l) 

(3) 

where e is the BER for the S-R link. 

IE. BER Analysis for A-MRC 

From (3), we observe that the BER has to be 
computed separately for the case of correct and 
incorrect decision at the relay. 



A. Correct Decision at Relay 

The probability of error, given a correct decision 
at the relay, can be expressed as 



P e n =P(X + AY < 0|x, = l,x r = 1) 
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F x (-Ay)p Y (y) dy 



(4) 



To obtain the above, the statistics of X and Y are 
required. Since X and Y are conditionally Gaussian 
[3], their statistics are known for integer values 
of the Nakagami fading parameters m,- [3], [4]. 
Using this, the BER in (3) was obtained in [4]. 
For arbitrary m,, while the exact PDF of X and Y 
is known (6), [4], an approximate expression for 
the CDF is available only for high SNR (7), using 
the Loskot-Prony approximation [6]. Due to space 
constraints, the proof of (7) is omitted in this paper. 
(4) can be expressed as 
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F x (-Ay)p Y (y) dy 
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F x (Ay)p ¥ (-y)dy (5) 



Substituting F x , a s > 0 from (6) and (7) in the first 
integral in (5), we have 
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F x (-Ay)p ¥ (y) dy = 
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Z I A ~ N(aA,bA),b > 0, A ~ @{c,m\ (a u a 2 ,a 3 ) = (0.168,0.144,0.002), (b u b 2 ,b 3 ) = 
(0.876,0.525,0.603), e = &,k = ^,K.(-) is the modified Bessel function of the second kind [7] and 



y (•, •) is the lower incomplete gamma function [7]. 
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The above integral is of the form 

I m<n (a,/3,6) = J y m+n e a >K m (Py)K n (6y)dy, 

{m,n,j3,6}>0. (9) 

This integral does not appear to be tabulated and 
is difficult to obtain in closed form. However, from 
(5), (6), (7), it is evident that the integral appears in 
the final expression for the BER and we will use (9) 
repeatedly in the following to represent integrals of 
the form in (8). The second integral in (5) can now 
be expressed as 
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The first integral in (10) can be expressed using 
integration by parts as 
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resulting in 
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upon substituting for F Y , a r > 0 from (7). From [9, 
(6.619.3)] 
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Using (13) in (12), we obtain 
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Thus, from (4)-(14), we obtain 
P(X + AY < 0\x s = l,x r = 1) = 
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5. Incorrect Decision at Relay 

Given that an incorrect decision is made at the 
relay, the probability of error can be expressed as 

P eh i =Fv(X + AY <0\x s = \,x r = -1) 
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The first integral in (16), upon substitution from (7) 
for F Y , a r < 0 is 
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where the third integral in (17) is obtained after 
subsituting for p x , a s > 0 from (6). The second 
integral in (17) can be expressed using the approach 
in (11) as 
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Substituting F x ,(a s > 0), from (7) in (18), 
1 

T(m r 



x mr l e^y(m s , -) 



^— f y(ra r , -?-\px(x)dx = 1 
m r ) Jo \ te r ) 

m r )\^£ r / I Jo 
Jo 



r(m,) 



x A I ^ Jh„ ih„ + j 



x mr l e^dx 



(19) 



The first integral in (19) can be expressed using [10, 
p. 138, (7)] as 
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The second integral in (19) can be expressed using 
(13) as 
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The third integral in (17) can be expressed using (9) 
as, 
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With this, all integrals in (17) are evaluated. The 
second integral in (16) can be expressed as 
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From (9), the second integral in (23) is obtained as 



From (16)-(24), we obtain 
P(X + AY<0\x s = \,x r = -1) 



3 0 /• 

_ ^ 2a i, 



b n 



ti T ( m r){e?(b n + * t )) 

2c n :° ( i 



*\a) T(m s )y/2nF s 



x 



+ J„ 



a s 2b„ 2 \~i ~T 1 , 



+ 



T{m r ) \Ae r ) j-f Y(m s ) 



3 

I 



a n 



X 



r(2m s + m r )Y(m r ) 



2m,+m r 



+ 



T(m s + m r + ^) 
2m. s + m r , m. s + i _ - ^ - f VM»» + 1))' 
ms+mr + I '(-1--^ + fVM^I)), 
1 / 6, \ mr 1 r(m s + m r ) 



r(m r ) T(m s ) 



m 



1 + 



x 2 Fi l,m v + m r ;m s + 1; 



1 + 



(25) 



2a„ 

Z^f(m^ 



1 / J \ /)! 



X 



2c?* 



r(m. s ) ^2nb s 



^a] + 2£. v c v , 



xJ„ 



(24) 



Substituting (15) and (25) in (3), we obtain the final 
expression for the BER. 

IV. Results 

In Figure 2, the analytical and simulated BER are 
plotted with respect to the average SNR for the S- 
D link. For convenience, we have chosen E r = E s , 
i.e. the source and relay transmit with equal power. 
(15) and (25) are used to compute the analytical 
BER using (3) for two cases, A = 0.5 and A = 1. 
As we can see, there is an excellent match between 
the simulation and analytical results, validating the 
expressions derived in the paper. Note that the 
Nakagami fading parameters are not integers. 

Figure 3 provides some interesting insights into 
the diversity order for /1-MRC cooperation. Firstly, 



V. Conclusions and future work 
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Fig. 2. Comparison of the simulation and analytical results for 
m = 3.7, m s = 2.6, m r = 2.6. Both match perfectly. 
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Fig. 3. Analytical BER plots for A = 1. Slopes for the lower two 
curves almost identical at high SNR indicating a similar diversity 
order. 



we note that the middle and bottom curves in Figure 
3 have the same slope at high SNR, indicating the 
same diversity order. We note that m s + m r = 2.7 for 
the middle curve is exactly equal to m = 2.7 for the 
bottom curve. This validates the result in [11] where 
the diversity order was shown to be min (m, m s + m r ) 
when A = 1. Note that the top curve has a diversity 
order 2 < 2.7 and its slope is less compared to that 
of the other two curves, at high SNR. 



In this paper, we have obtained a close but 
approximate expression for the BER for the A- 
MRC-DF cooperative system. The final expression 
contains only one integral in terms of simple, well 
defined functions. Numerical results obtained us- 
ing this expression match exactly with the actual 
simulation results, indicating the usefulness of this 
work. A closed form expression for the integral 
is a work in progress. Based on the techniques 
employed in this paper, it should be possible to find 
similar expressions for the BER for the superior PL 
combiner, and will be addressed in future work. 
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